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This work employs Hall magnetohydrodynamic (MHD) simulations to study the X- 
lines formed during the reconnection of magnetic fields with differing strengths and 
orientations embedded in plasmas of differing densities. Although random initial 
perturbations trigger the growth of X-lines with many orientations, at late time a 
few robust X-lines sharing an orientation reasonably consistent with the direction 
that maximizes the outflow speed, as predicted by Swisdak and Drake [Geophys. 
Res. Lett., 34, L11106, (2007)], dominate the system. The existence of reconnection 
in the geometry examined here contradicts the suggestion of Sonnerup [J. Geophys. 
Res. 79, 1546 (1974)] that reconnection occurs in a plane normal to the equilibrium 
current. At late time the growth of the X-lines stagnates, leaving them shorter than 
the simulation domain. 
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According to the frozen-in theorem of MHD, two adjoining collisionless plasmas with 
different densities, temperatures, and magnetic fields cannot alter their magnetic topology, 
and hence transport across their common boundary is prohibited. Magnetic reconnection 
violates this constraint as, for example, when solar wind plasma penetrates the Earth's 
magnetosphere or fusion plasma escapes from a tokamak core during a disruption. The 
questions of whether and how reconnection takes place for arbitrary plasma conditions are 
important for these and other systems. 

Consider two plasmas threaded by magnetic fields Bi and B2 of arbitrary relative orien- 
tation and separated by a planar surface, the x — z plane in Fig. [TJ We define the y axis to 
be perpendicular to the discontinuity plane, the z axis to parallel the X-line, and the x axis 
to complete the right-handed triplet. Let 9 be the shear angle between the two fields and 
a be the unknown angle between Bi and the X-line. In the highly symmetric cases often 
considered in theory and simulations a can frequently be easily deduced (e.g., a = 90° for 
Bi = — B2,). For more general configurations, however, no obvious choice exists, nor is it 
even clear that a single X-line orientation will dominate the system. Sonnerup- argued that 
a is determined by requiring that the currents in the reconnection plane vanish or, equiv- 
alently, that the components of the fields parallel to the X-line (the guide fields) be equal. 
As a consequence, no reconnection occurs in this scenario when cos 6* > B1/B2 (assuming, 
as in Fig. [TJ B\ < B 2 ), since no component of the field changes sign across the discontinuity. 

Others^ have questioned this choice on both theoretical and observational grounds. For 
example, although the Sonnerup criterion implies that reconnection between plasmas with 
small shear angles occurs infrequently, in situ observations in the solar wind suggest the 
contrary to be true^~-. In fact, most reconnection events in the solar wind occur at shear 
angles < 90°-^. 

As an alternative, Swisdak and Drake- proposed that the X-line orients itself so as to 
maximize the speed of the Alfvenic outflow. The outflow speed for plasmas with reconnecting 
components B\ x and B 2x and mass densities p\ and p 2 is&2 



Writing this expression in terms of a and maximizing with respect to a for a fixed 9 deter- 
mines the X-line orientation. Since v ont always has a local maximum between a = and 
a = 9 reconnection occurs for any 9 7^ 0. An alternative suggestion holds that maximiz- 




(1) 



2 



(a) 



(b) 



_ 

O 





► 




^ * V 1 






X 






.V 



FIG. 1. Field line geometries related to the Sonnerup— hypothesis. The reconnecting plasmas 
with fields B\ and B2 occupy the spaces y > and y < 0. In both panels the X-line parallels z, 
reconnection occurs in the x — y plane, and the fields are oriented such that the components parallel 
to the X-line are equal. Sonnerup— proposed that reconnection occurs when the x-components of 
B\ and B2 are anti-parallel (a), and that it otherwise does not (b). 



ing a related quantity, the normalized reconnection rate, determines the X-line orientation 
(M. A. Shay, private communication). 

In this work we perform 2-D and 3-D two-fluid simulations of reconnection between 
asymmetric plasmas in order to explore the generic development of X-lines. We first use a 
2-D simulation to demonstrate magnetic reconnection in a system with small shear angle 
with the caveat that, since 2-D simulations artificially impose the orientation of the X- 
line, studying the full development of the system necessitates a 3-D domain. Hence, we 
also consider a 3-D simulation of the same system. In previous investigations of 3D Hall 
reconnection^— , the initial configuration of anti-parallel fields confined nascent X-lines to 
one plane between the two plasmas. Initially localized X-lines grew in the direction of the 
electron current and, in some cases, extended over almost the entire computational domain. 
For the more general situation considered here, X-lines in the linear stage of development 
grow on different planes, known as rational surfaces, and undergo more complex interactions. 
We find that X-lines of several different orientations are excited at early times, but eventually 
only a few modes dominate. Interestingly, and in contrast to previous investigations^ 1 ^, 
the X-lines' length stagnates at a finite value that is shorter than the simulation domain. 
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For our initial equilibrium we employ a double tearing mode configuration with magnetic 
field components 

B x {y) = tanh ( y±hd±\ _ tanh ( V- L v^\ _ J (2) 

and 

B M = -4= tanh ( i±hl*\ + -L tanh ( ^Zi) + 2V2, (3) 

where w$ = 0.5 is the initial width of the current sheet (the normalization is described later). 
The asymptotic fields have components (B x , B y , B z ) = (1,0, y/2) and (— 1,0, 2a/2). Total 
pressure is balanced using a non-uniform number density n given by 

" = K f °-t)' (4) 

where the temperature T = 1.0 is uniform, B is the magnitude of the magnetic field, P a = 5.5 
is a constant, and the factor of 2 in the denominator arises from our code's normalization. 
For this configuration p 1 = 4, p 2 = 1, B\ = y/3, B 2 = 3, and the shear angle 9 = 54.7°. This 
system represents the limiting case described by Sonnerup- since cos# = B1/B2. 

The numerical simulations use the Hall-MHD code F3D^. It explicitly advances the 
dynamical variables (magnetic field, mass density, and ion velocity) with the second-order 
trapezoidal leapfrog method^ in time and fourth-order finite differencing in space. Periodic 
boundary conditions are applied in all directions. Variables are measured in normalized 
units: lengths to the ion inertial length di = (mjC 2 /47rnoe 2 ) 1 / 2 , velocities to the Alfven 
speed ca = So/(47rmjno) 1 ' 2 , densities to an arbitrary value no, pressures to Pq = minoc A , 
magnetic fields to an arbitrary field strength B , temperatures to rriiC 2 A and electric fields 
to Eq = caBq/c. Here c is the speed of light, and rrii and e are the mass and charge of the 
ions, respectively. 

The grid cells have a length of 0.2 on each side. No explicit viscosity or resistivity 
are applied, but a fourth-order diffusion coefficient of 10 -3 damps noise at the grid scale. 
The electron-to-ion mass ratio is m e /mi = 1/25. Since the electron inertial length d e = 
di-\Jm e /mi equals the cell size, the simulations do not describe the details of the electron 
dynamics. 

To determine whether reconnection can occur between the fields of Eqs. ([2]) and ([3]), we 
first perform a 2-D simulation. The computational domain has size L x x L y = 51.2 x 25.6 
and no variations are allowed in the z direction (i.e., d/dz = 0). We initiate reconnection 

4 



-40 -30 -20 -10 

X[dJ 

FIG. 2. (Color online) Out-of-plane current density J z for the lower sheet in the 2-D simulation at 
t = 156 with over-plotted magnetic field lines. The other sheet exhibits similar behavior. 

with random magnetic perturbations of amplitude 10 -3 £>o- The perturbations are generated 
in fc-space with maximum wavenumbers of k x = k y = 15. Setting the initial perturbation 
amplitude 100 times smaller produces similar final results. 

During the early stages of the simulation, strong out-of-plane currents develop, indicating 
the existence of magnetic reconnection. As the system evolves, multiple X-lines form, move 
along the x-axis, and merge, eventually leaving just one reconnection site. The out-of-plane 
current density of the lower current sheet after this merging is shown in Fig. [2j Following 
some initial fluctuations, caused by the interactions of multiple X-lines, the normalized 
reconnection rate stabilizes at a relatively steady value of ~ 0.02. 

This simulation demonstrates that reconnection can occur in a plane that includes equi- 
librium currents in a geometry where the model of Sonnerup- suggests it should not. How- 
ever, since the geometry of the computational domain determines the X-line orientation 
it did not establish whether an optimal orientation exists. Doing that requires a full 3D 
system, in which X-lines are free to develop in any direction. Our 3D simulation uses the 
same initial equilibrium as the 2D run, but applied on a computational domain of size 
L x x L y x L z = 51.2 x 51.2 x 409.6, and a grid scale in the out-of-plane direction of 0.2. 
Initial perturbations on the magnetic field in the z direction have a maximum wavenumber 
of k z = 5. 

In the linear theory of the tearing mode in periodic systems, reconnection can only occur 
at discrete locations, called rational surfaces, where k • Bo = (k is the wavenumber of 
the linear mode and Bq is the equilibrium field). Due to the periodicity of the domain, the 



wavenumbers of the linear instability must take the form k x = 2nm/L x and k z = 2irn/L z , 
where m and n are integers. This establishes rational surfaces at the locations satisfying 

B x (y) nL x 

B z {y) mL z U 
Since L z /L x = 8 several linear modes with n^O can grow in the current layer. Once they 
exit the linear regime interactions between different modes allow them to no longer respect 
the rational surfaces 

To reiterate, the model of Sonnerup- predicts that reconnection will not occur in this 
system. In contrast, Swisdak and Drake- predict that reconnection will occur with the 
X-line at an angle a given by the root of the equation 
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lying between and 9. Numerically solving for the parameters of the present system - 
B 2 /B 1 = v / 3,p 2 /pi = 1/4, and 6 = 54.7° — yields a = 34.3°. 

The early evolution of the 3-D simulation mirrors that of the 2-D run, in that both current 
sheets develop multiple X-lines separated by bulges in which reconnected flux accumulates. 
In Fig. [3] we present cuts at t = 84 (relatively early in the simulation) of J z in the x-y 
plane for four values of z separated by L z /4. For the 2-D simulation of Fig. [2] the X-lines 
necessarily parallel the z-axis, but that is not the case in Fig. [3j 

While in the 2D case topological constraints make the identification of X-lines straightfor- 
ward, in three dimensions the situation is more complicated, particularly when, as is the case 
here, no 3-D nulls exist in the system^. In this work we take an empirical approach to iden- 
tifying reconnection sites by examining isosurfaces of the current density J xz = \fj% + J z - 
(We use J xz in order to avoid favoring a particular X-line orientation in the x — z plane, 
although in practice we find that J z ^> J x .) Fig. [4] shows two values of the isosurface level 
at two different times. The regions of strong current, which map the X-lines, form extended 
structures in the x-z plane at various angles with respect to the z axis. At t = 84 (panels 
(a) and (b)) multiple X-lines are present, but by t = 201 (panels (c) and (d)) only a few 
distinct orientations dominate the system. Although the predominant X-line orientation at 
t = 201 is horizontal, there are some interesting features that appear to be aligned with the 
asymptotic magnetic fields. Although many of the weaker examples (for instance, the thin 

6 



51.2 

102.4 



X 



FIG. 3. (Color online) Slices of J z in the x — y plane of the 3-D domain at t = 84. The upper 
(lower) current sheet is described by bright (dark) colors, representing positive (negative) currents. 



structures at z ~ 40 and x ~ 20 in panel (c)) are field- aligned currents not directly associ- 
ated with reconnection, the strongest instances (e.g., the structure at —170 < z < —120 and 
— 10 < x < 10 in panel (d)) correspond to X-lines. Reconnection of the initial asymptotic 
fields cannot be occurring at such sites and, in fact, cuts through these features (not shown) 
reveal that the local reconnecting fields differ significantly from the initial asymptotic values. 

To evaluate the orientation of the X-lines quantitatively, we project the isosurfaces on 
the x — z plane and detect the edges with the Canny method^, a standard image processing 
tool that finds edges by looking for local maxima of the gradient. The gradient is calculated 
using the derivative of a Gaussian filter. The method uses two thresholds to detect strong 



and weak edges and includes the weak edges in the output only if they are connected to 
strong edges. In Fig. |5l the edges of the isosurface projection of Fig. HU are shown in black. 
Due to imperfections in either the image data or the edge detector, there may be missing 
points on the desired curve. 

The grouping of the extracted edge features to determine the X-line orientation is done 
with the Hough transform^. For each image pixel and its neighborhood, the Hough trans- 
form algorithm determines whether an edge exists at that pixel. The pixels lying along the 
highest values of parametric lines represent potential lines in the input image. Small gaps are 
automatically filled, and the lines are identified while a threshold is applied so that only lines 
longer than that value are considered. The red lines in Fig. [5] are the identified lines of the 
Hough transform, and clearly map the X-line. By simply averaging the various orientations 
identified in Fig. |5j we find that at t = 201 the X-line is oriented at an angle of ~ 0.6 ±8.5° 
with respect to the z-axis, which corresponds to a S i m = arctan(l/-y/2) + <p m 35.9° ± 8.5°. 

We have demonstrated that magnetic reconnection can occur between magnetic fields of 
small shear angle, in particular in configurations in which there is no anti-parallel component 
of the magnetic field in the plane perpendicular to the equilibrium current. This result 
disagrees with the model of Sonnerup-. In the initial stages of the simulation, we see 
the growth of X-lines with multiple orientations. By late time, however, one direction 
predominates. The orientation of this X-line, a S i m = 35.9° ±8.5°, agrees with the prediction 
of Swisdak and Drake-, olsd = 34.3°. The existence of many short X-lines with differing 
orientations early in the simulation demonstrates that L z does not play a limiting role and 
hence that the size of the computational domain probably does not affect our results. 

The value of olsd predicted by Swisdak and Drake- is the angle that maximizes the 
outflow speed from the X-line when the reconnecting fields have their asymptotic initial 
values. The strong diagonal features in Fig. HU discussed above, on the other hand, are 
due to the reconnection of significantly perturbed fields. (The density contrast, pz/pi, 
remains essentially constant.) Applying the criterion of Eq. [6] to the perturbed fields yields 
orientations roughly consistent with those observed in the simulation. It is unclear, however, 
if there is some overarching reason why the local reconnecting fields would be reconfigured 
in such a way as to generate X-lines that parallel the original asymptotic fields, or if the 
alignment in this case is purely coincidental. 

It has been suggested (M. A. Shay, private communication) that maximizing the normal- 
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FIG. 4. (Color online) Top view of the isosurfaces of J xz in the upper current layer at t = 84 
(panels (a) and (b)) and t = 201 (panels (c) and (d)). Panels (a) and (c) are at an isosurface level 
of 1.3, panels (b) and (d) at 1.7. The axes have been shifted to put the prominent features near 
the center. Solid lines denote the asymptotic magnetic fields and the dashed line the expected 
orientation of the X-line according to Swisdak and Drake—. Note the different axis scales. 



ized reconnection rate, and not the outflow speed from the X-line, determines the orientation. 
If the aspect ratio of the diffusion region (R, assumed to be < 1) remains independent of 
the upstream properties of the plasma (which has not been established in 3-D simulations) 
Cassak and Shay- argue that the normalized rate E varies as 

E ~2R^±( BlxB2x V (7) 



c V B\x + B 2x _ 

In Fig. E]the solid and dashed lines trace the dependence of v out and E ona for the parameters 
of our simulations; the vertical red line gives a s \ m . The dashed line peaks at an angle 
o>e = 31.7°, which is less than ccsd and slightly farther away from the value measured in the 

9 




-200 -150 -100 -50 50 100 150 200 

Z[d] 

FIG. 5. (Color online) Edges of the projection of the strong current density at t = 201 on the x — z 
plane and detected by the Canny method are shown in black. The red lines are the result of the 
Hough transform. 

simulation. However, given the relatively broad peaks generated by both criteria and the 
uncertainties associated with determining a S im, we cannot reliably discriminate between the 
two. The similarities between the quantities being maximized means that doing so probably 
requires extreme choices of parameters (e.g., B 2 /Bi 3> 1) that are difficult to simulate. 

As can be seen in Fig. HI the lengths of the X-lines barely change between t = 84 and 
t = 201; in fact most growth occurs early in the simulation, before significant magnetic 
flux has reconnected. We find that the growth of a given X-line is usually throttled by the 
interaction of its current with islands of reconnected flux from other X-lines at different 
rational surfaces. This effect will not be present in anti-parallel reconnection (where all of 
the X-lines are confined to a single plane) and may explain why our result conflicts with 
the finding of Huba and Rudakov— that X-lines continually grow in the current direction. 
Shay et alM did see stagnation of the X-line length for some initial current sheet widths, 
although not for the value used in this work (wo = 0.5). 

Our results suggest that reconnection can occur in any system where the adjoining fields 
are not parallel and in which other processes do not suppress reconnection (e.g., diamagnetic 
drifts^). The relatively broad peak of v ont in Fig. |6] may mean that, for a given set of 
asymptotic conditions, X-lines do not take on a single orientation but instead exhibit a 
distribution of orientations. Further 3-D simulations are needed to test this hypothesis. 

We are not aware of any other model that does a better job of predicting a. We suggest 
that on an encounter with reconnection events in which highly asymmetric conditions exist, 
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FIG. 6. (Color online) The magnitude of v ou t (Equation [TJ and the normalized reconnection rate 
of Cassak and Shay— (see Equation [7]) versus a. The ordinal units are arbitrary and have been 
suppressed. The vertical red line shows a s { m from our three-dimensional simulation. 



or while numerically reconstructing such an event, the Swisdak and Drake- criterion can 
cautiously be applied to determine the orientation of the reconnection X-line, as has already 
been done by Phan, Gosling, and Davis- and Teh and Sonnerup-. 
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